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Abstract 



It is shown that double inflation (two minimally coupled massive scalar fields) 
can produce correlated adiabatic and isocurvature primordial perturbations. 
Depending on the two relevant parameters of the model, the contributions 
to the primordial perturbations are computed, with special emphasis on the 
correlation, which can be quantitatively represented by a correlation spec- 
trum. Finally the primordial spectra are evolved numerically to obtain the 
CMBR anisotropy multipole expectation values. It turns out that the exis- 
tence of mixing and correlation can alter very significantly the temperature 
fluctuation predictions. 
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£j : I. INTRODUCTION 
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In our present picture, cosmological fluctuations today are seen as the combination of an 
initial spectrum, which can be computed within the framework of high energy models (like 
inflation), with the subsequent processes occuring at lower energy where the physics is be- 
lieved to be better understood (up to some unknowns such as the amount and nature of dark 
matter, the mass of neutrinos...). In the near future, we expect new and precise informa- 
tion on cosmological fluctuations with the planned measurements of the Cosmic Microwave 
Background Radiation (CMBR) anisotropics by the satellites MAP and PLANCK @. It 
has been emphasized in the recent years that the precision of these measurements could in 
principle allow us to determine with a high precision the cosmological parameters || . These 
studies however all assume very simple initial perturbations, typically Gaussian adiabatic 
perturbations with a power-law spectrum. However, reality could turn out to be more sub- 
tle. This then would have the drawback to complicate the determination of the cosmological 
parameters but could open the fascinating perspective to gain precious information on the 
primordial universe. At present, at a time when data are still unprecise, it is essential to 
identify broad categories of early universe models and to determine their specificities as far 
as observable quantities are concerned, with the purpose to be able to discriminate between 
these various classes of models when detailed data will become available. 



Ultimately, inflation must be related to a high energy physics model. Today there are 
many viable models but a generic feature of these models is that they contain generally 
many scalar fields. A property of inflation with several scalar fields is that it can generate, 
in addition to the ubiquitous adiabatic perturbations, isocurvature perturbations. In this 
respect, it is important to consider the possible role of primordial isocurvature perturbations. 
Isocurvature perturbations are perturbations in the relative density ratio between various 
species in the early universe, in contrast with the more standard adiabatic (or isentropic) 
perturbations which are perturbations in the total energy density with fixed particle number 
ratios. Primordial isocurvature perturbations are, most of the time, ignored in inflationary 
models. The main reason for this is that they are less universal than adiabatic perturbations 
because, on one hand, they can be produced only in multiple inflationary models [§J, and, on 
the other hand, they can survive until the present epoch only if at least one of the inflaton 
fields remains decoupled from ordinary matter during the whole history of the universe. 
However, not only the existence of isocurvature perturbations is allowed in principle, but 
candidates for inflatons with the required above conditions even exist in many theoretical 
models (dilatons, axions). 

What has already been established is that a pure isocurvature scale-invariant spectrum 
must be rejected because it predicts on large scales too large temperature anisotropies with 
respect to density fluctuations ||. But other possibilities can be envisaged. Several have 
been investigated in the literature: tilted isocurvature perturbations ||, combination of 
isocurvature and adiabatic perturbations j7|. In the latter case, only combinations of inde- 
pendent isocurvature and adiabatic perturbations were considered. The aim of this paper is 
to investigate the possibility of correlated mixtures of isocurvature and adiabatic perturba- 
tions. 

To illustrate this, the simplest model of multi-field inflation is considered here: double 
inflation ||, namely a model with two massive scalar fields without self-interaction or mu- 
tual interaction (other than gravitational). The production of fluctuations in this model 
has been studied in great detail by Polarski and Starobinsky and, in the present work, their 
notation and formalism will be followed closely. They were interested essentially in adiabatic 
perturbations [[| (see also |IOj for a numerical analysis) but also considered isocurvature per- 
turbations fll ]. However, they did not investigate the range of parameters where this model 



has the striking property to produce correlated isocurvature and adiabatic perturbations. 
By this, we mean the cases where both isocurvature and adiabatic perturbations receive 
significant contributions of at least one of the scalar fields, in contrast to the uncorrelated 
case where one of the scalar fields feeds essentially the adiabatic perturbations while the 
second one is at the origin of the isocurvature perturbations. 

The plan of this paper is the following. In section 2, the model of double inflation 
will be presented. Section 3 will be devoted to the analysis of adiabatic and isocurvature 
perturbations: their definition, how they are obtained from the inflation perturbations, 
the conditions to obtain correlated mixtures. Section 4 considers formally the definition 
of spectra for the perturbations as well as the notion of correlation. In section 5, the 
predictions for the CMBR anisotropies and matter power spectrum are given for the models 
with correlated primordial perturbations. 



II. DOUBLE INFLATION 

As mentioned in the introduction, inflation needs at least two scalar fields to produce 
isocurvature perturbations. That is why we investigate the simplest model of inflation with 
two scalar fields: they are non-interacting, massive, minimally coupled scalar fields. The 
Lagrangian corresponding to this model is 

( ^R 1 1 1 1 

C= 16^G~ 2 W " < ~ 2 m < < ~ 2 d ^ hdfl<Ph " 2 ml ^ (1) 

where the subscripts I and h designate respectively the light and heavy scalar fields (and 
thus rrih > mi). ^R is the scalar spacetime curvature and G Newton's constant. 

A. The background equations 

In a spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) spacetime, with met- 
ric ds 2 = —dt 2 + a 2 (t)<ix 2 , the equations of motion read 

3H 2 = 4ttG(0 2 + 4L + m 2 4> 2 + m 2 h 4> 2 h ) (2) 

4>i + 3H<j)i + m 2 (j)i = 0, 
'<j) h + 3H<j) h + m 2 J> h = 0. (3) 

Following it is convenient, during the phase when both scalar fields are slow-rolling (i.e. 
when (f) 2 and 4> 2 h can be neglected in (|J), (pi and (fit in (§)), to write the evolution of the two 
scalar fields in the following parametric form 



2 ^sinS, <Pi = ^cos6 (4) 

where 

s = — ln(o/a e ) (5) 

is the number of e-folds between a given instant and the end of inflation. This form (|]) is 
a consequence of the approximate relation d{(j>\ + 4> 2 )/ds = —d((f) 2 h + <ft 2 )/(Hdt) ~ [2-kG)^ 1 
resulting from the (slow-roll) equations of motion. The angular variable 9 can then be related 
to the parameter s by the expression 

2 

(sin^)^ 53 ! 

* = s »~ ST (6) 

(cos^)^ 2 -! 

where R is the ratio of the masses of the two scalar fields 

R = ^. (7) 

mi 

The equation @ was obtained by integrating the relation giving d9/d(\n s) as a function of 
6, which can be established by use of the slow-roll approximation of the equations of motion 



(0)-(H) (see || for the details of the calculations). As noticed in 0, equation (|j) can be 
approximated, when 9 3> -R -1 , by the simple formula 

(8) 



cos 2 6 



This behaviour corresponds to the period when inflation is dominated by the heavy scalar 
field (this approximation is valid as long as s > s and s — s ^> s /R 2 ). This period ends 
when 6 ~ -R -1 , s ~ s and is followed (possibly after a dust-like transition period) by a 
phase of inflation dominated by the light scalar field. 

It follows from @-(§j) that the Hubble parameter can be expressed in the form 

H 2 (s) ~ -m 2 s[l + (R 2 - 1) sin 2 6(s)], (9) 

o 

where the function 9(s) is obtained by inverting (g) (0 < 6 < n/2). As inflation proceeds, 
s decreases and 9 goes to smaller and smaller values, which implies a decreasing Hubble 
parameter during inflation. 

It will be convenient to define sh as the number of e-folds before the end of inflation 
when the scale corresponding to our Hubble radius today crossed out the Hubble radius 
during inflation. The value of sh depends on the temperature after the reheating (see e.g. 
jl^l ) but roughly sh — 60. To make definite calculations, we shall take throughout this 
work the value 

s H = 60. (10) 

Note that the class of models considered here depends on three free parameters: the two 
masses m\ and rrih, or alternatively mi and R, and the parameter s - in particular, the choice 
of this last parameter so relatively to s# will determine the specific phase of double inflation, 
'heavy' dominated, intermediate or 'light' dominated, during which the perturbations on 
scales of cosmological relevance were produced. 

B. Perturbations 

After having determined the evolution of the background quantities, let us turn now 
to the evolution of the linear perturbations. We shall restrict our analysis to the so-called 
scalar perturbations (in the terminology of Bardeen |I3||). We thus consider a spacetime 
linearly perturbed about the flat FLRW spacetime of the previous subsection, endowed with 
the metric 

ds 2 = -(1 + 2$)dt 2 + a 2 (t)(l - 2m)5 ij dx i dx j . (11) 

Although this metric is not the most general a priori, it turns out that any perturbed metric 
(of the scalar type) can be transformed into a metric of this form by a suitable coordinate 
transformation. This choice corresponds to the so-called longitudinal gauge. In addition to 
the geometrical perturbations $ and \1/, one must also consider the matter perturbations, 
which will simply be during inflation the perturbations of the scalar fields, respectively 8(f>h 
and S(pi, with respect to their homogeneous values. 



Before writing down the equations of motion for the perturbations, it is convenient to 
use a Fourier decomposition and to define the Fourier modes of any perturbed quantity / 
by the relation 

A= /t^W^/M- (12) 



(2vr) 3 / 2 

The equations of motion for the perturbations are derived from the perturbed Einstein equa- 
tions and from the Klein-Gordon equations of the scalar fields. They lead to the following 
four equations (see e.g. fll4"f ) 

$ = #, (13) 

$ + #$ = AnG (<j> h 6<j> h + frSfr) , (14) 



>h 



+ 3H6J> h +[^ + m 2 h \S(f> h = 40 h $ - 2m 2 h <f> h $, (15) 



5b + 3H6<h + I - + m 2 J ^ = 4<&$ - 2mf&$, (16) 

where the subscript k is here implicit, as it will be throughout this paper. 

In the slow-rolling approximation and for superhorizon modes, i.e. k <C aH, these 
equations can be solved (see 0) and the dominant solutions read 

$ ^ —F7T + 2C 3 V * 2 J '- h 2 h J 2 l , (17) 

Hi „,Ci nr Hml4> 2 h 5<p h „(h Hm 2 <f) 2 

where Ci(k) and C^k) are time-independent constants of integration and are fixed by the 
initial conditions. As usual in inflation, perturbations are assumed to be initially (i.e. before 
crossing out the Hubble radius) in their vacuum quantum state. Perturbations outside the 
Hubble radius are then obtained by amplification of the vacuum quantum fluctuations due to 
the gravitational interaction. The two scalar fields being independent, one simply duplicates 
the results of single scalar field inflation (see e.g. J14l). Consequently, 5(ph and 54>i can be 
written, for wavelengths crossing out the Hubble radius, as 

5<p h = J^=e h (k), 64*= -%=«(#, (19) 

where e^ and ei are classical Gaussian random fields with (e*(k)) = 0, (ej(k)e*(k')) = 
5ij5(k — k'), for i,j = l,h, and Hk is the Hubble parameter when the mode crosses the 
Hubble radius, i.e., when k = 2naH. Neglecting the evolution of the Hubble parameter 
with respect to that of the scale factor, the number of e-folds Sk corresponding to the 
instant when the mode of wavenumber k crossed out the Hubble radius, is given simply by 

k~k H e SH - s \ (20) 

where ku is the wavenumber corresponding to the present Hubble scale. In the present work, 
our interest will focus on the scales of cosmological relevance, typically k/kn — 0.1 — 2000. 
This means that the range of e-folds that will interest us is 52 < Sk ^ 63. 



III. PRIMORDIAL PERTURBATIONS 

The analysis of the solutions for the perturbations during inflation obtained in the pre- 
vious section will now enable us to determine the "initial" (but post-inflationary) conditions 
for the perturbations in the radiation era taking place after inflation and reheating. 

A. Initial conditions in the radiation era 

At some past instant deep in the radiation era, we shall consider four species of particles. 
Two species will be relativistic: photons and neutrinos; two species will be non-relativistic: 
baryons and cold dark matter. Their respective energy density contrasts will be denoted S 7 , 
8„, S b and S c (5 A = Sp A /p A ). 

At this point, it is useful to define precisely the notion of adiabatic and isocurvature 
perturbations. Isocurvature pertubations are defined by the condition that there is no 
perturbation of the energy density in the total comoving gauge (denoted by the subscript 
(c)), i.e. 

E^W = 0, (21) 

A 

but that there are perturbations in the ratios of species particle numbers, i.e. 

6 {c \n A /n B )^0 (22) 

in general. By contrast, adiabatic (or isentropic) perturbations are defined by the prescrip- 
tion that the particle number ratios between various species is fixed, i.e. 

S^ c \n A /n B )=0, (23) 

whereas the total energy density perturbation can fluctuate, i.e. 

£* (C W0. (24) 

A 

It is clear from the above definitions that, if one considers N species, there will be in general 
one adiabatic mode and N — 1 isocurvature modes. Here, it will be assumed that the light 
scalar field <pi decays into ordinary particles, i.e. gives birth to the photons, neutrinos and 
baryons, while the dark matter particles are associated exclusively with the heavy scalar 
field <ph- Note that part of the dark matter could be also produced by the light scalar field 
but we shall ignore this possiblity here for simplicity. As a consequence, the particle number 
ratios between the three 'ordinary' species will be frozen, i.e. 



5^ku L _ 6®n v _ S^n b 
n 7 n u rib 



(25) 



and only one isocurvature mode will exist, which can be conveniently represented by the 
quantity 



5 {c) n c _ ^ c >n 7 = (c) _ 3 



5 - — - — = ^ ~ A c) - ( 26 ) 



one can 



Going back to the longitudinal gauge, and following Ma and Bertschinger [jig 
write the initial conditions deep in the radiation era for modes outside the Hubble radius in 
the form 

S 1 = -2$, 

3 
5 C = S + ~8 1 , 

6 7 = e u = 6 b = e c = l(k 2 v)$, 

i 



2 



°v = i=(kr]) 2 $, 
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* = (l + -i^) $, (27) 

with R v = p„/(p-/ + p v ) and where 6 stands for the divergence of the fluid three- velocity, 
o v for the shear stress of neutrinos (in the rest of this paper, the contribution of neutrinos 
in all analytical calculations will be ignored for simplicity but it will be taken into account 
in the numerical calculations) and rj is the conformal time defined by dr] = dt/a(t). All the 
information about the initial conditions is thus contained in the two k-dependent quantities 
$ and S, which are time-independent during the radiation era (for k <^C aH). The next sub- 
sections will be devoted to make the link between these two quantities and the perturbations 
during inflation. 

B. Adiabatic initial perturbations 

The evolution of $ is given by ([T7|) only during the phase when both scalar fields are 
slow rolling. As soon as the heavy scalar field <fih ends its slow-rolling phase, the second term 
on the right hand side of ( |I7| ) will die out, and $ can be given during all the subsequent 
evolution of the universe in the simple form 

$ = Ci(l- — fait^dA. (28) 

It can be checked that, during inflation, 1 — -/ 'a(i')(lt' ~ —H/H 2 , which ensures that the 
coefficient C\ in the above formula is the same as in (|i~7|). 

It is thus essential to express C\ in terms of the perturbations of the two scalar fields, 
and therefore in terms of e\ and e^, in order to be able to determine the amplitude of the 
perturbations after the end of inflation. Combining the two equations in flT8"| ) and using 
(|19l), as well as the slow-roll approximation of the background equations of motion (0)-([|), 
the expression for the coefficient C\ during inflation is found to be 



d(k) ~ -4ttG-=|= [M(k) + <j> h e h (k)] . (29) 




As it is clear from this formula, C\(k) is a stochastic variable, whose properties can be 
determined from the stochastic properties of e\ and e^. 

During the radiation era, the relation between the coefficient C\ and the gravitational 
potential is simply, using once more (p8|) , 

$ = |d(k). (30) 

Therefore, the gravitational potential during the radiation era for modes larger than the 
Hubble scale is given by the expression Q 

$ ~ -^—k-^^Hu [sin 9 k e h (k) + cos fc e,(k)] , (31) 



where H k is given as a function of s^ in (H) and Sfc is given as a function of fc in (p0|). The 
hat in the above equations (and in all subsequent equations) indicates that the value of the 
corresponding quantity is taken deep in the radiation era when the wavelength of the Fourier 
mode is larger than the Hubble radius. 

C. Isocurvature initial perturbations 

As explained in subsection A, the isocurvature perturbations in the present model are 
due to variations in the relative proportions of cold dark matter, generated by the heavy 
scalar field, with respect to the three other main species (photons, baryons, neutrinos), 
all generated by the light scalar field. Moreover, during the radiation era, the isocurvature 
perturbation S rigorously defined by (P6|) is essentially the comoving cold dark matter density 
contrast, S ~ Sp, so that what is needed to obtain the primordial isocurvature spectrum 
is simply to compute the cold dark matter density contrast in terms of the scalar field 
perturbations during inflation. This task was carried out in ||11|| . Only the main points will 
be summarized here. 

Let us first give the comoving energy density perturbation associated with the heavy 
scalar field: 

Sp^ = <j) h 5<p h + m\<\> h b<\> h + 3HJ> h 6<p h - <Pl$. (32) 

Matching the inflationary phase when <ph is slow-rolling to the the inflationary phase when 



!>/j is oscillating and then to the post-reheating radiation dominated phase, one finds [J I 



5i c) * -\mlC (33) 

The coefficient C^ can then be obtained, during inflation, by substracting the two equa- 
tions in ( |18D and then using the (slow- roll) background equations of motion and (|T9|). In- 
serting the result in (|33|), the density contrast of the cold dark matter (associated with the 
heavy scalar field), for modes larger than the Hubble radius, is found to be given, during 
the radiation era, by the expression 

4 C) ~ ^H k (Ve^k) - ^k^(k)) , (34) 



where the value of the scalar fields is taken at Hubble radius crossing. This can be reex- 
pressed, using (|) and (]?]), in the form 



S ~ 5 { h c) ~ 2V^Gk- 3 / 2 sf /2 H k 



e h R 2 

-ei 



sin 9k cos 61. 



(35) 



Note that the isocurvature perturbations have the same power-law dependence as the adia- 
batic perturbations multiplied by a weakly k-dependent expression which is different from 
the analogous expression in flSTf). 

D. Conditions for the existence of correlated adiabatic and isocurvature 

perturbations 

As shown above, the quantities describing the primordial adiabatic and isocurvature 
perturbations are in general linear combinations of the independent stochastic quantities e; 
and eh and are thus expected to be correlated. It is now necessary to examine the actual 
value of the corresponding coefficients. For adiabatic perturbations, i.e. in equation (|3T|). 
the light contribution is dominant for tan 9 < 1 whereas the heavy contribution is dominant 
for tan 9 > 1. For isocurvature perturbations, i.e. in equation (|35|), the light contribution 
dominates for tan# > R 2 whereas the heavy contribution is predominant in the opposite 
case. Assuming R 2 3> 1, one can thus divide the space of parameters for double inflation 
into three regions: 



1. Region tan# >> 1 

The adiabatic perturbations are dominated by the heavy scalar field while the isocurva- 
ture perturbations are dominated by the light scalar field. The two types of perturbations 
will thus appear independent. Moreover, except for 9 very close to n/2, the isocurvature 
amplitude will be suppressed with respect to the adiabatic amplitude by a factor s^. In this 
parameter region, one recovers the standard results of a pure adiabatic spectrum due to a 
single scalar field, here the heavy scalar field. 

2. Region tan « R~ 2 

In this region, essentially only the light scalar field contributes to the adiabatic perturba- 
tions while the isocurvature perturbations are dominated by the heavy scalar field. The two 
contributions are therefore independent and the isocurvature amplitude can be very high 
with respect to the adiabatic one if 9 is sufficiently small. 

3. Region R- 2 < tan 9 < 1 

This is the most interesting region. Here, both the adiabatic and isocurvature perturba- 
tions are essentially feeded by the fluctuations of the light scalar field (even if their amplitude 
depends on the background value of the two scalar fields). This means that the adiabatic 



and isocurvature perturbations are strongly correlated in this region. If one considers the 
relative magnitude of these light scalar field contributions, one sees that the isocurvature 
contribution can compensate the s^ 1 suppression (with respect to the adiabatic perturba- 
tions) by a suitable factor R 2 . Note also that in the upper part of this parameter region, 
i.e. for tan# ^$ 1, the heavy and light contributions in the adiabatic perturbations will 
be of similar order while the heavy contribution in the isocurvature perturbations can be 
ignored. In contrast, in the lower part of the region, i.e. 9 ~ R~ 2 , the heavy contribution 
in the adiabatic perturbations is negligible whereas the light and heavy contributions in 
the isocurvature perturbations are of similar weight. This is illustrated on Figure 1, which 
displays the relative behaviour of the four contributions as a function of the angle 9. 

The expressions for the adiabatic and isocurvature primordial perturbations, (|3l|) and 
(|35|) can be simplified further when one assumes that these perturbations are produced 
during some specific phases of inflation. For instance, if all scales of cosmological relevance 
are produced during the period of inflation dominated by the heavy scalar field, then the 
approximate relation (|8]) relating the angle 9 to the number of e-folds applies, which enables 
us to simplify the Hubble parameter expression, given in equation (0), into 



H{s) ~ J -rmVW^lV^^- (36) 

The various contributions to adiabatic and isocurvature perturbations then reduce to the 
form 

k 3/2 ® h ~ -^—mi\/R 2 - 1 fe - s ) , k 3/2 ^ t ~ ^miVR 2 - lyjs (s k - s ), 

(37) 
and 



k*<*S h ^ 2 ^m l VW=l, k*"* - -^ mi E 2 ^S (38) 

3 3 y Sq 



where the indices h and I refer to the corresponding coefficients of e^ and e/ in (|31|) and 
(B5I). Let us briefly comment these results when one varies the free parameters of the 
model, mi, R and So (but remaining in the domain of validity of the above approximate 
expressions). Considering the variations with respect to the first two parameters, one can 
notice that all the contributions are proportional to the term mi\/R 2 — 1, except Si which 
contains an additional R 2 dependence. This means, ignoring for the moment the (weak) 
scale dependence, that the relative amplitudes of three of the contributions are fixed, the 
relative amplitude of Si being adjustable by the mass ratio R. Once R is fixed, the overall 
amplitude of the perturbations can be fixed by the scale mi. Concerning now the variation 
of the contributions with the cosmological scale, Sh is scale- invariant, while the three other 
are weakly scale dependent: $; and Si have the same dependence, whereas $^ has a stronger 
dependence. 

Another limiting case corresponds to 9 <C -R -1 , which occurs during the period of inflation 
dominated by the light scalar field. In this case, one has s ~ sq9 2 ^ r and the Hubble 
parameter is approximately given by 

10 




#0)~y~7W^. (39) 

As a consequence, the 'heavy' and 'light' contributions are approximated by 



k 3/2 $ h ~ ^—mis k 9 k , k 3/2 ^ L ~ -^m lSk (40) 

y y 



for the adiabatic perturbations and 



**»& = ^^Mji, ^5^-^p?m,fl 2 (41) 

for isocurvature perturbations. The 'heavy' adiabatic contribution is thus negligible and the 
perturbations due to the heavy scalar field are therefore essentially isocurvature. 

Note, to conclude this section, that in their work |J, Polarski and Starobinsky, concen- 
trated their attention on the intermediate case where the scales of cosmological relevance 
just correspond to the transition zone from the heavy scalar field driven inflation to the light 
scalar field driven inflation. As a consequence, their spectrum has a stronger variation in k 
than in the limiting cases considered above. Here, the emphasis is put on contributions to 
the isocurvature perturbations. With another choice of parameters, one can also produce a 
huge temperature anisotropy dipole due to isocurvature perturbations on scales larger than 
the present Hubble radius |L6| . 



IV. SPECTRA AND CORRELATION OF ADIABATIC AND ISOCURVATURE 

PERTURBATIONS 

A. General definitions 

It is usually assumed in cosmology that the perturbations can be described by (homoge- 
neous and isotropic) gaussian random fields. In the specific model under consideration here, 
where the perturbations are created during an inflationary phase, this is true by construc- 
tion. What is new here is that isocurvature and adiabatic perturbations are not assumed to 
be independent. Indeed, as shown in the previous section, in the case of double inflation, the 
two kinds of perturbations are correlated, at least for some region of the parameter space. 
It will thus be our purpose to define statistical quantities that can describe random fields 
which are, a priori, correlated. Let us first recall, for any homogeneous and isotropic random 
field /, the standard definition (up to a normalization factor) of its power spectrum by the 
expression 

(hf:,} = 27r 2 k- 3 V f (k)6(k-k'). (42) 

In addition to this definition, it will be useful to define a covariance spectrum between two 
random fields / and g by the following expression 

TleifaZ,) = 27r 2 k~ 3 C Lg (k)5(k - k'). (43) 
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In order to estimate the degree of correlation between two quantities, it is convenient to also 
define the correlation spectrum Cf^(k) by normalizing Cf t9 (k): 

^ {k) = /JnflnV (44) 



Schwartz inequality implies, as usual, that — 1 < Cf )9 (k) < 1. The correlation (anticorrela- 
tion) will be stronger as one is closer to 1 or — 1. 

B. Double inflation generated perturbations 

Let us now specialize the above formulas to the case of perturbations generated by 
double inflation. By substituting the explicit expressions for the perturbations obtained in 
the previous section, namely (|31|) and ( j35|) , one finds 



n = 


^ HkSk 


itrum, 


S * S k 


i? 4 1 
cos 2 9 sin 2 9 



(45) 

for the initial adiabatic spectrum, 

on n? T r 4 i 1 

(46) 
cu»" o biii a 

for the initial isocurvature spectrum and 

AG 
C*,s = -^H 2 k (R 2 - 1) (47) 



for the covariance spectrum. Combining the three above spectra according to (|44{), one finds 
finally for the correlation spectrum the expression 

*' s 2(R*sm 2 6 + cos 2 ey/ 2 ' [ ' 

It is instructive to study the dependence of this correlation spectrum with respect to the 
parameters of the model. If one takes 9 fixed, one sees that the correlation will vanish for 
R = 1 and will then increase monotonously with increasing R approaching the asymptotic 
value cos#. If now one considers R as fixed and study the variations of the correlation with 
respect to 6, one recovers the conclusions of section 3D: the correlation vanishes when 8 
approaches zero or 7r/2; inbetween, one can see that the correlation reach a maximum for 
sin 2 ^ = (R 2 + l) -1 , with the value 

n2_ i 

The correlation spectrum C^ § for various choices of parameters has been plotted on Fig. 2, as 
a function of s&. One can, as before, distinguish between the two extreme cases. For models 
such that 9 ^> R~ l , corresponding to a 'heavy' inflationary phase, the various contributions 
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vary slowly with 9, as can be seen from Fig. 1. This means that the correlation will be 
almost constant. For models such that 9 <C -R -1 , corresponding to a 'light' inflationary 
phase, §h increases quickly with decreasing 9, i.e. with decreasing scales, which implies that 
the correlation will decrease with decreasing scales, i.e. smaller s^. The models with s < sh 
belong to the first category, while models with sq > sh correspond to the second. Finally, 
the models with so close to sh have an intermediate behaviour between the two extreme 
cases. They also have the strongest correlation. 

V. PREDICTIONS FOR THE CMBR AND DENSITY CONTRAST SPECTRUM 
A. Analytical predictions for long-wavelength perturbations 

1. Evolution of the perturbations 

In the case of perturbations whose wavelength is larger than the Hubble radius, the time 
evolution is particularly simple. For an initial isocurvature perturbation characterized by 
the initial amplitude S, the entropy perturbation S is unchanged as long as the perturbation 
is larger than the Hubble radius, whatever the evolution of the backgroung equation of state, 
i.e. 

S = S {k^aH). (50) 

However, the radiation-matter transition will generate a gravitational potential perturbation 
(see e.g. 0) 

& so = --S (k<^aH). (51) 



Of course, the initial adiabatic perturbation will also contribute to the gravitational potential 
perturbation: 

$ od = T$ (k < aH), (52) 

where T is a coefficient, close to 1, due to the evolution of the universe (if one ignores the 
anisotropic stress of the neutrinos, T = 9/10.) 

2. Large angular scale CMBR anisotropics. 
At large angular scales, the temperature anisotropies are essentially due to the sum of 



an intrinsic contribution and of a Sachs- Wolfe |17| contribution. Except for the dipole for 
which the Doppler terms are important, the Sachs- Wolfe contribution can be written (for a 
spatially flat background) 

/AT\ 1 

(^r) (e) = ^M; (53) 

V 1 J SW o 

Where e on the left hand side is a unit vector corresponding to the direction of observation 
and Xi s on the right hand side represents the intersection of the last scattering surface 
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with the light-ray of direction e. The intrinsic contribution is simply given, via the Stefan 
law, as the perturbation #p/4 at the time of last scattering. Since last scattering occured 
in the matter era, 6$ — 5^ c \ and therefore for an adiabatic perturbation (5$ = jS!^), 
(-^r). — \$mi which can be seen to be negligible (see below (|56|)) with respect to the 

Sachs- Wolfe contribution, whereas for an isocurvature perturbation (S ~ — |<5j^ during 
matter era), 

J- J int O 

To conclude, the temperature anisotropies will be given in general by 

*I = ir* - 2 -s (55) 

T 3 5 K ' 

on angular scales larger than the angle (of the order of the degree) corresponding to the size 
of the Hubble radius at the time of the last scattering. This equation enables us to estimate 
easily the normalization of the temperatures anisotropies for the low multipoles (see the 
definition fl62|)-(|63|)), essentially constrained by COBE measurements. Note that for mixed 
primordial perturbations with isocurvature and adiabatic contributions of the same order of 
magnitude, the low multipoles anisotropies can be significantly reduced by a compensation 
effect between the isocurvature perturbation and the adiabatic one. It turns out that this 
is the case for the light scalar field contribution in double inflation models with R ~ 5 (see 
Fig.l and the consequence on Fig. 3-5). 



3. Large scale structure 

Large scale structure is governed by the density contrast, or equivalently the gravitational 
potential perturbation $ since the latter quantity can be related to the (total) contrast 
density in the comoving gauge by the (generalized) Poisson equation, which reads (see e.g. 

0) ' 

= ~\^\ (56) 

For modes inside the Hubble radius the evolution becomes quite complicated and depends 
on the specific ingredients of the model. But what is relevant for our purpose is that this 
subhorizon evolution does not depend on the nature of the primordial perturbations. What 
matters is the total gravitational potential perturbation $ which can be written, in the 
matter era, as 

$ = <5> ad - ~S. (57) 

5 

Note that the influence of primordial isocurvature perturbations is smaller on the large scale 




density power spectrum (see (0)) than on large scale temperature anisotropies (see (|55"D) 
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4- Spectra 

Using the relation fl55p (with T = 1), the spectrum for the large scale temperature 
anisotropies can be expressed in terms of the primordial isocurvature and adiabatic spectra, 



=r V* + 25^ 15' 



V ^ = *n + ~^s - ~^s- (58) 



When only primordial adiabatic perturbations are present, the previous expression implies 

V'Jl = \vT, (59) 

T 6 

whereas for pure isocurvature perturbations, one finds 

V l Jl = 2V X J 2 . (60) 

T 

This is in agreement with the standard comment in the literature that isocurvature perturba- 
tions generate CMBR anisotropies six times bigger than equivalent adiabatic perturbations. 
This is the reason why isocurvature perturbations are in general rejected in comological 
models 0. However, when one takes into account both isocurvature and adiabatic pertur- 
bations, with the possibility of correlation, then the additional term due to correlation can 
change significantly these conclusions. Illustrations will be given in the next subsection. 
Similarly, the spectrum for the gravitational potential is given by 

?* = ?* + ^P S - \C^. (61) 

B. All-scale predictions 

After having considered long wavelength perturbations, whose advantage is one can es- 
timate analytically their observable amplitude and thus normalize easily the models, let us 
analyze now smaller scales, which require the use of numerical computation. 

1. CMBR anisotropies 

As it is customary, one decomposes the CMBR anisotropies on the basis of spherical 
harmonics: 



AT 



oo I 

E E a»m*im(M). (62) 

1=1 m=—l 



The predictions of a model are usually given in terms of the expectation values of the squared 
multipole coefficients 

Q = (h m | 2 >. (63) 
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In the present model, the temperature anisotropies will be the superposition of a contribution 
due to the heavy scalar field and of a contribution due to the light scalar field. These two 
contributions are independent, because the stochastic quantities e^ and e/ are independent, 
and as such the coefficients Ci can be decomposed 

Cl = ci l) + Ci h \ (64) 

where the upper index refers to the 'light' or 'heavy' nature of the perturbations. It is 
important to emphasize that only a decomposition of this type is allowed here. For example, 
a decomposition of the C\ as a sum of an isocurvature contribution and of an adiabatic 
contribution would be wrong here. In practice, the heavy and light contributions to the C% 
are computed independently, by using twice a Boltzmann code (developped in our group 
by A. Riazuelo, and used in flT8|). The first run takes as initial condition $/j and S/, and 



yields the coefficients C\ . Similarly, the second run computes the C\ using as initial 
conditions the corresponding quantities ^ and Si. The results for C\ and C\ , as well as 
their sum Ci, are plotted on Fig. 3-6 for four illustrative models (for all models the Hubble 
parameter and the baryon density correspond respectively to /iioo = 0.5, Qb = 0.052). For the 
first three models, the value R = 5 has been chosen because the isocurvature and adiabatic 
contributions of the light scalar field are then of similar amplitude, as is visible on Fig. 1, and 
the effects of mixing and correlation are particularly important. A consequence of the similar 
amplitude (with the same sign) of the two 'light' contributions is an important suppression 
of the light spectrum C\ for small /, as noticed already in the previous subsection, and as 
is visible on Fig 3-5. In contrast, one can check on Fig. 6 that this will not be the case for 
the R = 10 model, for which Si is dominant. 

The first two graphs have a roughly similar behaviour for the 'heavy' contribution. What 
distinguishes them is the 'light' contribution, which illustrates its high sensitivity on the 
relative amplitudes of 5/ and $;. A systematic investigation of the effects of mixed correlated 
primordial spectra, independently of the early universe model to produce them, on the 
temperature anisotropies will be given elsewhere ]HJ. For these two models, one notices an 
amplification (weak in the first case and strong in the second) of the main acoustic peak 
with respect to the standard (pure adiabatic and scale- invariant) model. In contrast, the 
third example shows a suppression of the main peak, which is due to a strong contribution 
Sh, which makes the 'heavy' spectrum look "isocurvature" and thus damps the main peak 
in the global spectrum. Finally, the last example is characteristic of the domination of the 
'light' spectrum, itself dominated by the isocurvature contribution (Si), which thus makes 
the global spectrum look "isocurvature". It is rather remarkable that modest variations 
of the two relevant parameters of the model, R and Sq (mi is useful simply for an overall 
normalization of the parameters), can lead to a large variety of temperature anisotropy 
spectra. 

C. Power Spectrum 

Another quantity which is extremely important for the confrontation of models with 
the observations is the (total) density power spectrum. In the literature, it is usually de- 
noted P(k) and its relation to the corresponding spectrum V S ( C ) (for the comoving density 
constrast) defined generically in (|42]) is 
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P{k) = 27r 2 k- s V 5(c) . (65) 

Using the Poisson equation (|5(]), it can be reexpressed in terms of the gravitational potential 
spectrum 

p{k) = iSw kn{k) - (66) 

As with the temperature anisotropics, the power spectrum for double inflation is obtain- 
able by computing independently the power spectrum for the heavy scalar field contribution, 
then that for the light scalar field contribution, and finally by adding the two results, 

P(k) = P l (k)+P h (k). (67) 

In contrast with the temperature anisotropies, the influence of the mixing and correlation of 
the primordial perturbations on the density spectrum is less spectacular because the shapes 
of the pure isocurvature and pure adiabatic density spectra are not extremely different. 
There is however a sensible difference: the pure isocurvature spectrum has, relatively to 
the large scales, less power on small scales than the pure adiabatic spectrum. To illustrate 
what happens for mixed and correlated primordial perturbations, Fig. 7 displays, for the 
model specified by the parameters R = 5 and s = 50, the total power spectrum , together 
with the two independent 'light' and 'heavy' contributions, as well as the standard adiabatic 
CDM power spectrum (normalized as before) for comparison. Note that the resulting spec- 
trum has, relatively to large scales, less power than the standard adiabatic power spectrum. 
However, it has globally more power than the standard spectrum for the same temperature 
anisotropy amplitude on small I. 



VI. CONCLUSIONS 

The main conclusion of this work is that it is possible, in the simplest model of multiple 
inflation, to obtain correlated isocurvature and adiabatic primordial perturbations. These 
perturbations slightly deviate from scale-invariance but their correlation can entail significant 
modifications with respect to standard single scalar field models. 

This class of models, both simple and rich, could provide an interesting field of experiment 
to investigate the feasibility to determine the cosmological parameters and the primordial 
perturbations from expected data. The question would be, assuming Nature has chosen 
this particular model, could we infer from the expected temperature anisotropy data the 
cosmological parameters and to which precision? More important, would it possible to 
discriminate between a single field inflation model and a multiple field model with correlated 
perturbations and what would be the price to pay on the precision of the cosmological 
parameters ? 

It was not the purpose of the present work to exhibit a model supposed to fit better the 
observations. However, one of these models, surprisingly, turns out to present two character- 
istics, which are at present favoured by observations: a power spectrum with modest power 
at small scales (comparatively to the standard CDM model) and a high peak on intermediate 
scales. It may be worth seeing how well this model does when confronted with the current 
observations. 
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Finally, it would be interesting to investigate the possiblity of correlated adiabatic and 
isocurvature perturbations within the framework of multiple inflation with interaction be- 
tween scalar fields and see how the main features presented here would be modified. 
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FIG. 1. Relative amplitude of the (frh contributions to the adiabatic perturbations, <I>/j (dashed 
line), to the isocurvature perturbations, Sh (dotted dashed line), and of the (pi contribution to the 
adiabatic perturbations, <3?/ (continuous line), and to the isocurvature perturbations, Si (dotted 
line). The last quantity, the only one which depends on R, has been plotted for R = 5 (upper 
dotted line) and R = 10 (lower dotted line). 
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FIG. 2. Correlation spectrum for various parameters. Continuous curves from bottom to top 
(on the left hand side of the figure) correspond respectively to (R = 5, so = 30), (R = 5, so = 50) 
and (R = 10, so = 50). The dashed curve corresponds to (R = 5, so = 60), the dotted dashed curve 
to (R = 5, sq = 70) and finally the dotted curve to (R = 5, So = 80). 
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FIG. 3. Temperature anisotropies for the double inflation scenario with R = 5, sq = 30. The 
total anisotropies (continuous line) are the sum of a contribution due to the heavy scalar field 
(dashed line) and of a contribution of the light scalar field (dotted line). To make the comparison, 
the anisotropies due to standard (adiabatic scale-invariant) perturbations are also plotted (dotted 
dashed line), using Cio for normalization. 
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FIG. 4. Temperature anisotropies for the double inflation scenario with R = 5, so = 50 (same 
conventions as in Fig. 3). 
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FIG. 5. Temperature anisotropics for the double inflation scenario with R = 5, sq = 80. 
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FIG. 6. Temperature anisotropies for the double inflation scenario with R = 10, sq = 50. 
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FIG. 7. Power spectrum for the double inflation scenario with R 



5, so = 50. The total power 
spectrum (continuous line) is the sum of a contribution due to the heavy scalar field (dashed line) 
and of a contribution of the light scalar field (dotted line). The standard (adiabatic scale-invariant) 
power spectrum is also plotted (dotted dashed line) with the same normalization for Cio- 
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